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Resumen 

En este trabajo se presenta un modelo matemático empirico para la función de Gibbs 
de exceso, g E =y/(p,T,x), de un sistema multicomponente. La dependencia con la composición se 
realiza a través de la denominada fracción activa, z¡, la cual a su vez, se relaciona con la frac- 
ción molar jc¡ de los componentes de una solución y de un parámetro k y] cuya determinación se 
indica en el trabajo. Se analiza Ia eficacia del modelo en base a su extensión, considerando 
varios casos y tres posibilidades de cálculo para el parámetro k X] , Io cual da lugar a un conjun- 
to de versiones del modelo para la correlación de datos cityas ventajas se discuten. El modelo 
propuesto para la función de Gibbs de exceso adopta la expresión genérica siguiente, 

g E (P ,7>)=z(*)[l-z(jc)]Xgi0 P,T)z■' siendo g¡(P.T)=g a +g l2 P 2 +g ii PT+^-+g i5 r 

i=0 1 

y permite su aplicación a un caso general de equilibrio líquido-vapor con variación de las tres 
variables fundamentales jc¡, p y T, o bien teniendo en cuenta valores experimentales de dos ejer - 
cicios importantes, isobárico e isotérmico, que también se analizan. La derivación matemática 
permite obtener otras propiedades de mezcla, e incluso llevar a cabo una correlación simultá- 
nea de varias de ellas. El modelo se ha aplicado a varios sistemas binarios para los que existen 
datos experimentales en la literatura en un amplio rango de p y T. Los resultados puede consi- 
derarse aceptables. 

Palabras clave: Propiedades termodinámicas, Función de Gibbs, Equilibrio de fases, Correla- 
ciones. 


Summary 

In this paper an empiricial mathematical model for excess Gibbs function, g=y/(p, T.x). 
of a multicomponent system, is presented. The dependence with the composition is carried out 
by using of the concept of active fraction, z¡, which is related with the molar fraction x\ of the 
components of a solution and with a parameter k\- } the determination of which is indicated in this 
paper. The eficiency of the model is analyzed considering several aspects such as its extension 
with several cases, and three possibilities to calculate the parameter k X] showing a set of possi- 
bilities for the correlation of data. All cases are analyzed and discussed. The model proposed 
for the excess Gibbs function adopts the following generic expression, 

siendo g i (P.T)=g n +g i2 P'-+g 0 PT+^+g i5 r- 

i=0 1 



and permits its application to a general case of vapor-liquid equilibria considering the varia- 
tion of three fundamental variables x lf p y T, or taking into account experimental values of two 
relevant problems, such as the isobaric and isothermic cases, which are also analyzed. The 
mathematical derivative permits geting other mixing properties and even the simultaneous cor- 
relation of several of them. The model has been applied to several several binary systems for 
which experimental data are available in the literature in a wide range of p and T. Results can 
be considered acceptable. 

Key words: Thermodynamic properties, Gibbs function, Phase equilibria, Correlation. 


1. Introducción 


Teniendo en cuenta que actualmente no existe una teoría general y satisfactoria 
que proporcione una base sólida para la estimación, o para la correlación de propiedades 
termodinámicas de sistemas multicomponentes, puede considerarse al trabajo empírico 
como una forma de acercamiento a lo primero, lo cual no significa que se abandonen los 
esfuerzos en esa dirección. Durante varios años se han realizado publicaciones donde se 
presentaron ciertas expresiones empíricas que pueden proporcionar de forma más dire- 
cta los valores de algunas cantidades termodinámicas de mezcla de sistemas binarios 
cuando se conocen las medidas experimentales de las mismas. En ellas se ha detallado 
la representación matemática de esas propiedades mediante funciones empíricas que no 
presentan relaciones conocidas con la física molecular o la mecánica estadística. 

Para contribuir a la modelización de soluciones, las dos propiedades de mayor 
interés, tanto bajo un punto de vista práctico como teórico, son las cantidades de exceso 
de la energía de Gibbs y la entalpía; la entropía y otras pueden calcularse a partir de 
estas. Un complemento a ellas sería también el análisis de los cambios de volumen de 
mezcla, cuya determinación es bastante precisa con los medios actuales a partir de la 
medida directa de la densidad. Por ello, en este trabajo centramos especialmente nuestra 
atención en la función de energía de Gibbs de exceso, ya que, a partir de ella pueden 
obtenerse las demás cantidades termodinámicas de las soluciones y facilitar su interpre- 
tación. En cualquier obra de Termodinámica Química pueden conseguirse las siguientes 
relaciones para la función de Gibbs G, que permiten conseguir, a una determinada com- 
posición, otras propiedades. 
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La utilización de estas relaciones con las correspondientes cantidades de exceso 
resulta inmediata. Estas expresiones demuestran el enorme interés que tendría establecer 
una única ecuación o modelo para la función de Gibbs de una solución homogénea, vá- 
lida para condiciones diferentes de presión y temperatura. Tan importante es su deter- 
minación experimental como su estimación teórica y el tratamiento matemático que se 
lleve a cabo con esos datos, ya que la inter/extrapolación, cuando no existen datos sufi- 
cientes, resultaría necesaria en muchos casos. 
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2. Representación de los datos de equilibrio. Diagrama de fases 

La aplicación de la regla de las fases a un sistema en equilibrio constituye un 
importante recurso para la construcción e interpretación de su diagrama de fases a través 
del conocimento del número de variables independientes de un determinado sistema. La 
regla de las fases, dictada por J. Willard Gibbs en 1875, establece precisamente el nú- 
mero de variables independientes, o grados de libertad, que deben fijarse de manera 
arbitraria en un sistema para definir su estado intensivo. La mencionada regla se expresa 
matemáticamente por: L=C+2-F, donde L es el número de variables independientes, C 
el número de componentes y F el número de fases del sistema en estudio. Conociendo el 
número de variables independientes pueden establecerse también las restricciones geo- 
métricas sobre los aspectos que se requieren para describir las fases coexistentes en un 
diagrama. Según dicha regla, el número máximo de variables independientes para un 
sistema de 2 componentes en equilibrio es de 3. Por ello, el comportamiento del equili- 
brio de ese sistema puede describirse completamente por medio de volúmenes, superfi- 
cies, líneas y puntos en un espacio tridimensional. 

Por ser el caso que trataremos aquí, de sistemas binarios en dos fases, la aplica- 
ción de la regla de Gibbs a un sistema de dos componentes A y B y dos fases (1) (líqui- 
do) y vapor (v) en equilibrio, da lugar a: 

Grados de libertad=L= X[p, T, x A , v A ] - Y[//L v ; ;/ l =// v j = 4 _ 2 = 2 (3) 

En esta regla de Gibbs, las llamadas variables independientes son principalmente 
los potenciales, p, Ty p x , porque la expresión anterior se basa precisamente en la ecua- 
ción de Gibbs-Duhem que se planteará en este estudio. En lugar de potenciales pueden 
emplearse cantidades molares para cada una de las fases y que están relacionadas con 
los potenciales mencionados, a pesar de ser intensivas, sin embargo dichas cantidades 
no son potenciales. Con la relación (3), la definición de potencial químico y la condi- 
ción de equilibrio planteada por la ecuación de Gibbs-Duhem resuelve el planteamiento 
matemático de un problema real. Sin embargo, es importante obtener respuestas numé- 
ricas al modelo así planteado. Para ello, es necesario relacionar el potencial químico con 
el resto de las variables fundamentales escritas en el primer sumando de la ecuación (3), 
es decir con la p, T y ,y¡. Este aspecto está fuera del ámbito de la Termodinámica en el 
sentido más estricto, pero que es un problema a resolver. E1 conocimento de otras fun- 
ciones menos abstractas que el potencial químico, como la fugacidad y la actividad no 
ayudan a encontrar la solución. En resumen, las reglas de Gibbs permiten conocer la 
ecuación temiodinámica del equilibrio y emplearla como relación de partida (igualdad 
de potenciales químicos para cada compuesto en las diferentes fases en que está presen- 
te), sin embargo, es necesario conocer la relación de los //¡ con el resto de variables, p, 
T, A*j. Incluso, bajo un punto de vista gráfico sería adecuado representar las variables 
potenciales para cada compuesto en cada fase y encontrar la línea de coexistencia. En 
este sentido, este trabajo supone una contribución inicial dando un primer paso para 
conseguir relaciones con esas pretensiones, conociendo de antemano que la Termodi- 
námica no es suficiente para resolver el problema planteado. Por ello, un planteamiento 
empírico puede ser de a>aida. 

A continuación se presenta un breve análisis de los aspectos prácticos del caso 
que nos ocupa para un equilibrio líquido-vapor de un sistema binario. La ecuación (3) 
nos da 2 grados de libertad y las dos variables independientes medidas usualmente en 
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los equilibrios son py T, sin embargo, es necesario medir una tercera variable para des- 
cribir completamente el sistema de dos componentes y dos fases. Para conseguir una 
descripción cualitativa y comparar los datos de equilibrio, la elección más conveniente 
de la tercera variable es, bien la composición de la fase líquida, o la de vapor. Es decir, 
del conjunto de las cuatro variables, p, T, x,y, que intervienen, se necesita conocer tres 
(una condición puede ser la de trabajo isobárico o isotérmico). A pesar de ello, es difícil 
predecir teóricamente la cuarta a partir de las tres medidas previamente. Por ello, gene- 
ralmente la composición de la otra fase también se mide. La medidas de las cuatro va- 
riables p, T, x, y , si bien implica una sobreespecifícación en términos de la regla de las 
fases es, no obstante, el procedimiento preferido para una presentación precisa de los 
datos. Además, esta sobreespecifícación permite claramente la verifícación de los datos 
experimentales con la ecuación de Gibbs-Duhem. Lo que parece claro, al menos expe- 
rimentalmente, y a ello se refíere la regla de las fases de Gibbs, es que, para un caso 
isobárico por ejemplo, al fíjar la presión y una de las variables, v, la dinámica del siste- 
ma en ebullición a p=Cte da lugar a una temperatura de equilibrio T y una composición 
de la fase vapor, también invariables. 

En cuanto a la representación del 
comportamiento de fases, las condicio- 
nes de equilibrio en un sistema de dos 
componentes y dos fases puede descri- 
birse completamente en un diagrama 
3D con ejes para p, T, x, y, ver Figura 1. 
Sin embargo, dicha presentación de los 
datos de equilibrio es tediosa de dibujar 
y difícil de interpretar. Por esta razón, 
para visualizar en mejor forma el equi- 
librio de fases, se prefíere dibujar las 
secciones obtenidas cuando se mantiene 
constante una de las variables, ver Fi- 
guras 2(a-d). Si elegimos como varia- 
bles a las composiciones de un com- 
puesto en las dos fases, las demás per- 
manecen fíjas, lo que implica que las 
fracciones molares de ambas bases de- 
pendan de p y de T, es decir, que tam- 
bién serán x{p,T) y y{p,T). Estas funcio- 



Figura 1. Diagrama 3D de fases líquido-vapor 
estándar oara un sistema de dos comnonentes 


nes se corresponden a dos superfícies en el espacio p-T-x, o p-T-y , dando lugar, por 
ejemplo, a un corte transversal en la Figura 1 a p=C te, ver Figura 2(a). E1 signifícado de 
esa gráfíca es importante por partida doble; por un lado, por que da relaciones T-x de las 
fases en equilibrio, y en segundo lugar porque la zona de líquido y de vapor divide al 
plano T-x en tres regiones, la superior de vapor donde las mezclas de vapor tienen la 
más baja energía de Gibbs, lo mismo para la zona inferior de la zona líquida, y la región 
entre liquido y vapor que corresponde a la zona de más baja energía de Gibbs realizada 
como una combinación de mezclas de vapor y líquido. 

En base a todo esto, un aparato experimental se diseña para medir los datos de 
cada una de las secciones transversales del diagrama de fases tridimensional. Más con- 
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cretamente, los aparatos se diseñan para hacer medidas en únicamente una de las sec- 
ciones. Esto no resulta problemático ya que únicamente un diagrama particular de fases 
se requiere para el diseño de una operación. Los datos que mayormente se presentan en 
la literatura son los isotérmicos p-x-y, sin embargo, las medidas de datos isobáricos T-.x- 
v son mayormente demandados que los isotérmicos ya que la mayor parte de las opera- 
ciones básicas que se realizan en el campo de la ingeniería química se llevan a efecto a 
presión constante. E1 diagrama a composición constante p-T , tiene un empleo más limi- 
tado. En la Figura 2 (a-d) se han representado los diagramas parciales de cada uno de 
los planos de la Figura 1, además de la que corresponde a vi versus .vi, para un sistema 
binario (l)+(2) estándar. 






Figura 2.- Secciones del digrama p-T-x-y. (a) diagrama T-x-y para tres presiones; (b) diagrama p-x- 
v para tres temperaturas; (c) diagrama p-T para tres composiciones; (d) diagrama y-x a una presión 

E1 caso que nos ocupa en este trabajo podría parecer algo más complejo, ya que, 
además de las cuatro variables comentadas, todas ellas medibles experimentalmente, en 
los estudios de equilibrio entre fases puede trabajarse con una dependencia de las mis- 
mas interrelacionadas con otras obtenidas indirectamente, como es la función de energía 
de Gibbs. 

3. La función de Gibbs para una solución multicomponente 

De manera general, la fúnción de Gibbs de un sistema material constituido por 
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varias especies depende de la temperatura, de la presión, y del número moles de cada 
una de las especies presentes. Matemáticamente, esta dependencia funcional de la pro- 
piedad G puede escribirse como G=G{p y T,n x ), indicando con n\ al número de moles de 
cada una de las especies o sustancias que intervienen en un sistema multicomponente, es 

decir, de n\, m, .Explícitamente para una mezcla de C componentes puede escribir- 

se: 


G=G( p,T,n { ,n 2 ,...,n c ) 

Si derivamos la (4) con relación a cada una de las variables 
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pero como G es una funcion homogenea de orden uno con respecto a la cantidad de los 
componentes, esto es, G = gn, entonces: 


d( gn)= 


S(gn) 

dT 


dT + 


r d(gnT 


dp 


dp+I, 


/ ci(gn)"' 


/ T,n 


dn ; 


dn, 


( 6 ) 


i y 


P>T,nj*¡ 


E1 último sumando se corresponde con la propiedad molar parcial, G ¡9 por lo que la 
ecuación (6) también puede escribirse: 
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dp+Yfi y dn\ 


(7) 


T,x 


ahora, esta última ecuación, teniendo en cuenta las relaciones mostradas en (1) y (2), 
puede también plantearse por la expresión equivalente: 

d(gn)=-nsdT+nvdp+YG\dn\ (8) 

i 

pero: d(gn)=ndg+gdn y, dn x =x x dn+ndx \, entonces, al sustituir en (8), 


dg+sdT-vdp~YG ¡ dx\ 


g~TG\X\ 


dn=0 


(9) 


pero como n y dn son cantidades arbitrarias e independientes, una de las soluciones de 
la ecuación (9) es que ambos sumandos sean iguales a cero. Entonces puede plantearse, 
por un lado, 

dg=-sdT+vdp+YG [ dx l , (10) 


y por otro que g=YG\X\ » de donde dg=Yfi x dx ¡ +Y*\dG x , 


( 11 ) 


quedando entonces la ecuación (10) como sigue: 

-sdT+vdp-Yx\dG= 0 (12) 

¡ 

que se conoce con el nombre de ecuación general de Gibbs-Duhem. Para el caso espe- 
cial de p y T constantes, la ecuación (12) queda reducida únicamente al tercero de los 
sumandos, resultando entonces una forma restringida de la ecuación de Gibbs-Duhem, 
que indica que las propiedades molares parciales están acopladas. Regresando de nuevo 
a la ecuación (10) vemos que la función “g” resulta tener como variables a la presión p, 
la temperatura T, y la composición. Para el trabajo con soluciones, los dos primeros 
sumandos son establecidos a composición constante y pueden escribirse en función de 


120 






























otras cantidades termodinámicas medibles, mientras que el último de ellos recoge las 
correspondientes propiedades molares parciales de la propiedad con derivadas en 
relación al número de moles, n x , lo cual no resulta del todo práctico. Por ello, con el fin 
de ser rigurosos con el planteamiento físico del problema, intentaremos poner dicho 
ténnino en función de la composición o fracción molar de cada componente. 

Veamos ahora la iníluencia del cambio de variable en el término que contiene a 

G,. 
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y como /7=Z/7i, entonces {dn¡dn\) pTnj =1 • Ahora sería conveniente poner la propiedad 

intensiva “g” en función de las composiciones o fracciones molares .v¡, de los compues- 
tos que integran la solución, para las que se verifica Xxj=l. Así, para una disolución de 
“C” componentes, a presión y temperatura constantes, la propiedad depende de “C-l” 
fracciones molares, por lo que, eligiendo una de ellas como variable dependiente del 
resto, puede ponerse x¡=l-Xxj, y por tanto, 

g=g(x\,x 2 , ...,X¡_ 1 ,-Y¡+ 1 , ...,x c ) (14) 

ydeaquí: í/g=Zhp- | dx k (15) 

g T.xj.u 

pero para conseguir la derivada de la propiedad con relación al número de moles, ver 
(13), se necesita recordar la siguiente relación: x k =n k /n, de donde surge la derivada 
•Tk 


¡ =-Üíl=-^ vk=l,2.i-l.i+1.C 

M J p j, nj „ " 2 » 

teniendo en cuenta este resultado, en la ecuación (15) la derivada se plantea como: 
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y al reemplazar este término en (13) se produce 
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esta expresión puede fínalmente ser introducida en la ecuación (10) 
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g-X-V k dx¡ =-sdT+vdp+g Z<¿v¡ -ZZ- v k 

8x k J i , k U-v k 


(16) 


(17) 


(18) 


dx¡ (19) 


P,T,xj^u 


pero, de la relación entre la fracciones molares, se obtiene que S¿y¡= 0, por lo que la 
ecuación anterior se reduce fínalmente a: 
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Ecuación que contiene claramente las derivadas correspondientes a las variables pre- 
sión, temperatura y fracciones molares de los compuestos que integran el sistema. Por 
tanto, de manera simplificada y bajo un punto de vista matemático, esta ecuación puede 
asimilarse a la derivada de una función multivariable del tipo 
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T, X\ 9 X2, .)=0 (21) 

que ofrece una forma más adecuada para conseguir un modelo sencillo que pueda em- 
plearse en el tratamiento de datos y para la programación de los cálculos que la de parti- 
da, la ecuación (4). Las variables que aparecen en la posible ecuación empírica (21) son 
intensivas, independientes de la cantidad total de sustancia presente, aunque las x\ no 
son todas independientes entre sí, si lo son de T y de p. 

Una relación como la (21) también resulta al aplicar a la ecuación (4) una de las 
propiedades de las funciones homogéneas, al dividirla por el número total de moles n de 
la disolución, es decir: 

Grplj^ ,ii2 . >K l =G ( p j' Xí ' X2 . >Xc )=g{ p j iX] ' X2 . Xc ) (22) 


La diferencial de la ecuación (21) da lugar a una expresión cuya analogía matemática 
con la (20) resulta clara, que aun cuando alcanza el significado físico/termodinámico 
que allí se plantea, permite obtener una ecuación empírica adecuada para modelizar el 
tratamiento de datos como se ha comentado. La derivada de (22) será: 
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pJ,Xj*k.c 


Este objetivo se desarrollará en las siguientes secciones. Apuntar que todo este 
planteamiento será aplicable a la correspondiente magnitud de exceso, g E , en cuyo aná- 
lisis se concretan los estudios de equilibrio entre fases, por lo que, a partir de ahora, 
centraremos esta presentación en el desarrollo de dicha cantidad. 


4. Una expresión para C e (jc) 


Para implementar ahora un modelo de soluciones se debe realizar un ejercicio 
para defínir unas relaciones, y además con cierto sentido físico, de la función de energia 
de Gibbs de exceso con las variables de composición, presión y temperatura. Dichas 
variables han sido introducidas en trabajos anteriores [1-3] para situaciones particulares 
de diversas propiedades. La expresión más genérica del modelo tendrá la forma 



/>=2 


iii; 


Z a \{\ 2 ....\ p z \ x z \ 2 

i p <=CR*(n,p) 


(24) 


Así, se estableció que una cantidad de exceso viene definida como la suma de las con- 
tribuciones de los efectos de todas las posibles interacciones p-arias de las n fracciones 
efectivas, donde p= 2, 3, 4, etc. Q representa el máximo orden de las interacciones mole- 
culares que se consideren, mientras que los a x ¿ 2 -ip son los coefícientes particulares de 


influencia de la combinación p-aria en la cantidad de exceso que se considere; incluso, 
en algún caso uno o varios de esos coeficientes pueden ser nulos. Para la fracción activa 
z¡ se considera una relación con la fracción molar de la manera genérica siguiente: 

k " X ' (¡ = 2,..., rí) (25) 


Z: =~ 


- v i +Z^i- v í 


siendo ahora k\\ un coefíciente a determinar y al que dedicaremos un tratamiento sepa- 
rado, ver Anexos A.l y A.2. Cuando se desarrolla la ecuación (24) puede llegarse a una 
expresión recurrente para el caso de n sustancias, que puede escribirse genéricamente de 
la siguiente forma: 
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( 26 ) 


a/„ l o= 

il¡2 ■-in-1 €C(n.«-l) 

donde A/^'/q 2 "'~ ln_l) representa a cada frocción o contribución parciai de la cantidad 

de exceso debida a la interacción de (a 7- 1) componentes hasta un máximo de Q molécu- 
las elegidas entre las especies ii, L, ..., i„_i. Cada una de las contribuciones parciales a la 
función de exceso, puede a su vez reducirse a las A/^_ 1Q correspondientes, y así hasta 

llegar únicamente a las interacciones binarias. E1 término representa esencial- 

mente al producto entre z\zi ... z„ y un polinomio cuyo grado depende de Q. Para el caso 
de un sistema temario donde se consideren interacciones hasta de cuarto orden, la co- 
rrespondiente ecuación (26) se concreta en: 

A / 3 E 4 = A/ 2 e, 4 '' 2) + M, E( 4 '- 3) + + - 1 - 2-3 (Co +c,r, + C 2 z 2 ) (27) 

Aquí, los términos M E 4 ' _3 ' ) se refieren a la contribución parcial de la cantidad de exceso 

considerada del sistema binario i-j al cómputo global de la energía de mezcla del siste- 
ma temario con interacciones de cuarto orden como máximo. 

Como el conjunto de ecuaciones (24-27) puede aplicarse a cualquiera de las can- 
tidades de exceso, por la naturaleza de este trabajo la dirigiremos a la función de Gibbs 
de exceso g h . Para conseguir lo que se pretende, en lugar de utilizar una de las ecuacio- 
nes generalizadas, emplearemos la particular de un sistema binario simple, donde se 
cumple que xi+.y^-I y tambiénzi+Z 2 =l. De esta manera, la dependencia de la función de 
Gibbs de exceso con la composición, partiendo de la (24), y después de organizar térmi- 
nos y simplificar, queda reducida a: 

Ai- z iÍ a o +a i z í+ a 2 z ¡] ( 28 ) 

donde los coeficientes A¡ son una combinación lineal de los coeficientes particulares de 
influencia , ¡ 1Í2 ...¡ /7 . y, al menos aquí, independientes de la temperatura; la dependencia 

de la composición x¡ está clara a través de las z¡ [ecuación (25)]. 

5. Una expresión para G e (a% I , p) 

La cuestión que se plantea ahora es obtener un modelo único para la función de 
energía de Gibbs de exceso de un sistema multicomponente, esto es, lograr una expre- 
sión que contenga las tres variables principales, con una forma explícita del tipo 
g E =g E {P’T>x j). Por ello, en trabajos anteriores se plantearon las relaciones parciales 
que acogieran, en primer lugar, una dependencia únicamente de la composición .v [1,2] 
[ver ecuaciones (28) y (25)]. Luego, en otro [3] la correspondiente dependecia de la 
composición y la temperatura a través de: 

^+Wl-z^z' donde: A, =¿4^' + AJ* (29) 

/v 1 /=0 

y por último, de la composición y la presión [4]: 

¿M=z(l-z)¿S,z‘ siendo: B,=fb,p> (30) 

K1 

Con todo lo anterior, la función que se busca será una expresión en fiinción de la 
concentración, similar a la planteada para los dos casos estudiados, es decir, una fiinción 
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polinómica de la fracción activaz, pero donde ahora los coeficientes deberán contener a 
las variables intensivas de temperatura y presión. Replanteando en este instante una 
extensión de lo anterior para la función de Gibbs de exceso tenemos: 
g E (p,T,x) = z t z 2 [g 0 ( p, T) + g t (p, T)z { (x) + g 2 (p, T)z¡ (x)] (31) 

pudiendo ponerse a los gipj) como combinación de los coeficientes A¡ de (29) y B\ de 
(30), de la forma que se expresa a continuación y que contendrá a las expresiones par- 
ciales planteadas para las mencionadas ecuaciones. 


g¡ ( P’ T ')=g a +g a P 2 +g a pT+^jr+gj 2 


(32) 


En esta sección del artículo es importante matizar una importante cuestión sobre 
el modelo cuyo desarrollo se presenta. Tal como se ha expresado en las ecuaciones a las 
que hemos hecho referencia, en artículos anteriores se ha planteado el desarrollo de los 
modelos particulares (isobáricos o isotérmicos) para la función adimensional de Gibbs 
de exceso, g iRT , mientras que a partir de la ecuación (31), con la que se pretende esta- 
blecer un modelo más general que acoja a las tres variables p, T , x, se ha definido una 
expresión para la función de Gibbs de exceso en su forma intensiva y no-adimensional. 
Trabajar con una función adimensional se planteó en un principio porque en el cálculo 
de cantidades tennodinámicas donde se utilizan los ELV, es costumbre emplear datos 
de la función de Gibbs en su forma adimensional por su relación directa con los coefi- 
cientes de actividad. Las derivadas de la relación (32) con respecto a cada una de las 
variables p y T son: 


dg\(p.T) 

dT 


p,x 


=gnP-jY+ 2 g i5 T y 


( p,T ) 


dp 


=2pg\i+g\iT 


JT,x 


(33) 


siendo ambas ecuaciones del tipo que (29) y (30). Finalmente la siguiente es una posible 

forma desarrollada para la ecuación que buscamos: 

2 

g ,: ( pT,x)=z(\-z)Y J g i z'=z t z 2 [g <> +g t z+g 1 z 2 y 
¡=0 


=z,z n 


go\+go2P 2 +goiPT + ~- + gosT 2 +(g\\+g\ 2 P 2 +g\3PT + ^Z7 + g\sT 2 Z + 
i y V í ) 


g 2\ + g 22 P 2 + g 23 PT + -^" + £ 25 ^ 
donde la expresión rigurosa para z\ será: z¡ =—— (i = 2 ,..., n) 

Xx+2L k ÁP> T ) x \ 


(34) 


(35) 


Aún cuando la ecuación (34) parece tener un elevado número de parámetros, en la prác- 
tica dicho número es reducido, como se verá más adelante. Es importante conocer si 
esta expresión es consistente con todo lo planteado en artículos anteriores, siguiendo un 
formalismo matemático-termodinámico y más concretamente lo establecido para las 
relaciones (1) y (2). En primer lugar, es preciso determinar las derivadas de la función 
de Gibbs de exceso con respecto a la temperatura y con respecto a la presión y conse- 
guir así relaciones para otras propiedades, o sea que: 


h E =g E -T 


% 


E X 


dT 


-g V ~Tz\Z 2 Y, 


ydT 


dzY dk(p,T) '] 

dk A dT J 


(\-2z)£ gi z ¡ +z(l-z)£ig l z i - 


(36) 
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si para simplificar denominamos 
Y=(\-2z)Í gi z' + ;(l-z)Íig,z'-' 


( 37 ) 


1=0 


expresión que puede generalizarse como: K=^(i+lXg¡ -g,_, )z' siendo para nuestra 

► la expresiói 

1'=.- 0---)I 


aplicación g-i=g 3 = 0 . Introduciendo la expresión de (37) en (36 

\dT ) ^ {dT 


h^gt-Tz^^z'-T 


z'-T 


dz 


\dT j 


Y 


otras cantidades son: 




dg±\ 

dT 


(dY 


KdT 


-TY\ 




(38) 


(39) 


incluso, el término entre corchetes puede plantearse de una fornta más genérica, de la 
siguiente forma: 


análogamente, para la variación con p. 


V E = 

( E X 

dg h 

=z,r, y 

°gi ^ 

.-J*! 

Í«T 


{ d p j 

T.X ' 

l°P J 

T \dk) 

T,x v y 

l d P J. 


(*-2z)X^i zi +*0-*)£tei- i_ 


que podrá ponerse, considerando la ecuación (37), de la forma: 


= Z , Z :I 


( Z \ 

dg¡ 

\ d p j 


z + 


T.x 


dz_ 

dp 


y, de igual fonna, la función de Helmoltz será: 


(41) 

(42) 

(43) 

(44) 

(45) 


resultando en las ecuaciones (38) y (42) polinomios similares a los origínales para su 
empleo en el tratamiento de los datos de volúmenes de exceso [1] y de entalpías [2]. 

Por tanto, utilizaremos a partir de aquí el modelo básico (34) para su aplicación 
concreta a ciertos sistemas binarios para los que la bibliografia presenta datos de varias 
de las pro^iedades, en particular nos referimos a valores de ELV (isobáricos e isotérmi- 
cos), de h E , y de v E , también obtenidos a varias T y p. 


=Z \ = 2 X 


¡=o 


g\~P 


'dgA 

-« Jdz_J 

{ dk \ 

V dp J_ 

f 

% 

UJ. 


{\-2z)^g r .' + z(\-z)Y\gr-'- 

¡=0 7=0 


o bien de manera más simple: 



g\~P 

dg \ j 

i f d =) 

- -P\ — 

i=0 


dp J. 

\dpj 


6. Aplicación del modelo a sistemas binarios 

E1 modelo establecido en (34) puede ser aplicado a sistemas para lo que se co- 
nozca su comportamiento de fases en un amplio rango de temperaturas y de presiones, 
que puedan llegar a superponerse, o al menos que ofrezcan un cierto acercamiento en p 
y/o en T. A pesar del gran número de investigadores que actualmente trabajan en el 
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campo del equilibrio entre fases, no se encuentran demasiados sistemas para los que se 
conozcan valores experimentales de varias de las propiedades que se ponen en juego en 
el modelo propuesto. Una exhaustiva búsqueda bibliográfica ha indicado que dichos 
datos existen para algunos conjuntos o “grupos/familias” de las mezclas alcanol+alcano, 
dialquilcarbonato+alcano, ester+alcanol, y aromático+alcano. 

En la Tabla 1 se presenta un resumen de varias versiones del modelo propuesto, 
utilizando como forma base o estándar la presentada en (34). Se dan diferentes casos 
según se modifique, por un lado, el grado del polinomio en función de la variable com- 
posición, o de la fracción activa z, y por otro, el orden del polinomio de cada coeficiente 
de “z” como ftmción de las variables temperatura y presión. De esta forma se plantean 
tres grados diferentes para el polinomio y su vez, cada uno de ellos con cuatro o cinco 
coeficientes, lo que da lugar a los seis casos que se han calificado de “A” a “F” en la 
mencionada Tabla 1. Además, otra posibilidad que debe considerarse es la dependencia 
o no del parámetro k\¡ con p y T, hecho que también se analiza. 


Tabla 1. Relación de diferentes versiones del modelo, denominadas como A-F, en función del grado del 
polinomio, en z, y de la extensión de los coeficientes g x (p, T). 


modelo, gradoz 

ecuación 

coeficientesg¡(p,7) 

A, n=l 

g E =Z\Zl[gQ+g\Z\] 

g\( P.T)=g\\+gi2P 2 +gi3pT+g»/T 

B, n=l 

g E =Z\Z2[gO+g\Z\] 

g\(p.T )=g\ i +g\2P 2 +g\2PT+g\A/T+g\s/T 2 

C, n=2 

g E =Z\Z 2 [gO+g\Z\+g2Z\] 

g¡ < P.T )=g¡ 1 +g¡ 2 P 2 +gúpT +g¡4 /T 

D, n=2 

g E=z \Z2[gO+g\ z \+g2Z\] 

g\( P.T )=g\\+g\2P 2 +gi3pT+gi4/T+g\s/T 2 

E, n=3 

g V =2\Z 2 [gO +g\=\ +glP\+ g}-\ ] 

g\( P.T )=g\ | +g¡2P 2+ g\3P T+ 8\4/ T 

F, n=3 

g E ==\Z2 [gO +gl-1 +g2-\ +gl=\ ] 

g\( pT )=g\\ +g\2p 2 +g\$pT +g\4¡T+g\$/T2 


En cualquiera de los casos se utiliza una función objetivo ( F.O .) que engloba las 
diferentes cantidades o propiedades termodinámicas empleadas en la correlación de 
datos. La función objetivo se plantea como: 


n p N 


F 0.= 2 2 




i,exp / i,cal 


i,exp y 


P=y(iso-T),y(iso-p),h E ,v E ,Cp ; 


(46) 


siendo N el número de puntos, y n p el número de propiedades consideradas en la corre- 
lación. La Figura 3 recoge un esquema sobre el procedimiento global a seguir en la apli- 
cación del modelo y que detallaremos en los siguientes comentarios. 

Para llevar a cabo una aplicación práctica se distinguirán dos grandes casos. E1 
primero, utilizando el modelo de manera particular para, bien un proceso isobárico o un 
proceso isotérmico pero teniendo en cuenta las distintas versiones del modelo según la 
Tabla 1. Se establece para el parámetro k\¡ por ser, al menos en un principio, el que pro- 
duce una más óptima correlación, el que corresponde a los /r¡j-libre (nos referimos a k w 
para la correlación de los v E , a k^ para la correlación de las h ¥ \ y a k g para la correlación 
de las g E ). 
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Figura 3. Diagrama de flujo correspondiente a las vias/procedimientos que se plantean en la 
correlación de datos o valores experimentales con el modelo de la ecuación (34) y sus deri\ adas 


Con ello, el tratamiento de datos realizado en este primer escalón del proceso, 
dejando cada uno de esos valores libre en el procedimiento de ajuste, produce los resul- 
tados de la Tabla 2 para el conjunto de sistemas que allí se reflejan. Esta consideración 
del parámetro /ry-libre, utilizada para los seis modelos que se plantean en la Tabla 1 (A a 
F) y las desviaciones estándar obtenidas, dan lugar al conjunto de gráfícas, una para 
cada propiedad, de la Figura 3. En color ‘ k gris ,, se han sombreado las columnas corres- 
pondientes al modelo que da lugar a la desviación mínima de las propiedades conside- 
radas en conjunto. Ello produce una valoración cualitativa global de dicha aplicación, y 
es que la versión del modelo calificada como “D” y como “F v ocasiona las menores 
desviaciones estándar en casi todos los casos (en algunos puede incluirse al “C' y al 
’ k ET), si bien la primera de las versiones con un grado inferior del polinomio en “z \ por 
lo que, con el fin de evitar la sobreparametrización, consideraremos más idónea a efec- 
tos prácticos la versión Indicar que las desviaciones encontradas para las c l p no 

fiieron demasiado significativas en ninguna de las versiones establecidas, observándose 
en la Figura 4 una cierta uniformidad en los resultados de todos los modelos para cada 
uno de sistemas. Las mayores desviaciones para los sistemas 8 y 9 no son significativas. 
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Tabla 2. Error estándar, s , de los ajustes conjuntos de datos de ELV en condiciones isotérmicas e isobáricas, 
de entalpías de exceso, ^/(J-moÍ' 1 ), volúmenes de exceso, 10 9 v E /(m J -mor'), y capacidades térmicas de 
exceso, c* /(J-mol^K' 1 ) (cuando existen) mediante diferentes versiones del modelo (33) y califícados según 


la Tabla 1. Los parámetros k e , k h , k v , k c se obtienen en el procedimiento de ajuste. 


Sistema binario 

Error 



Modelo 



estándar 

A 

B 

C 

D 

E 

F 

Alcohol + Alcano 


1 

Etanol+r 2 Hexano 

S(g E )’ 

0.030 

0.031 

0.021 

0.018 

0.018 

0.017 


Ref.[5], Ref.[6], Ref.[7], Ref.[8] 

s(y-S 

2.769 

2.772 

1.760 

1.464 

1.241 

1.245 



S(/)" 

0.027 

0.027 

0.021 

0.016 

0.016 

0.015 



*0’¡)" 

1.689 

1.687 

1.166 

1.067 

1.058 

1.035 



s(h E ) 

98 

88 

86 

66 

82 

53 



5(v b ) 

33 

34 

17 

17 

17 

17 

2 

jqEtanol+r^Heptano 

í(g E )’ 

0.026 

0.028 

0.023 

0.024 

0.022 

0.023 


Ref.[5], Ref.[6], Ref.[7[, Ref.[8[ 

s(Yi)‘ 

0.412 

0.419 

0.301 

0.266 

0.224 

0.247 



S(/)" 

0.020 

0.020 

0.014 

0.010 

0.010 

0.009 




2.032 

2.030 

1.706 

1.548 

1.544 

1.562 



s(h k ) 

119 

109 

102 

75 

84 

52 



s(v B ) 

49 

49 

15 

15 

15 

16 

3 

x\ Etanol+Jc 2 Octano 

s(/)' 

0.047 

0.044 

0.048 

0.045 

0.047 

0.041 


Ref.[5], Ref.[6], Ref.[10] 

s(7i)' 

0.679 

0.666 

0.595 

0.525 

0.570 

0.462 



s(/)" 

0.015 

0.015 

0.011 

0.011 

0.011 

0.011 




0.439 

0.442 

0.388 

0.381 

0.372 

0.365 



s(h b ) 

127 

115 

112 

73 

105 

52 



s(v h ) 

47 

48 

13 

12 

13 

12 

4 

X[Etanol+x 2 Nonano 

s(g E )** 

0.029 

0.027 

0.031 

0.029 

0.029 

0.029 


Ref.[5] 

s(VÍ)" 

0.497 

0.514 

0.457 

0.423 

0.440 

0.422 



s(h E ) 

116 

113 

98 

68 

89 

54 



s(v h ) 

58 

58 

12 

13 

12 

12 

5 

X\ Btanol+XiDecano 

S(/)" 

0.073 

0.071 

0.071 

0.068 

0.070 

0.073 


Ref.[5] 

s(7i)" 

0.422 

0.449 

0.384 

0.418 

0.358 

0.313 



s(h E ) 

140 

130 

106 

66 

98 

80 



s(v h ) 

67 

65 

19 

20 

16 

15 

6 

X\ Propan-1 -ol+x 2 Heptano 

s(/)' 

0.019 

0.019 

0.013 

0.019 

0.019 

0.020 


Ref.[ll],Ref.[12], Ref.[13], 

s(Vi) 

0.332 

0.321 

0.156 

0.218 

0.218 

0.280 


Ref.[14] 

s(gT 

0.027 

0.027 

0.027 

0.030 

0.024 

0.033 



s(Vi)" 

0.188 

0.180 

0.232 

0.198 

0.255 

0.207 



s(h E ) 

126 

120 

102 

95 

79 

36 



s(v) 

20 

20 

7 

7 

7 

8 

Dialquilcarbonato + Alcano 


7 

X]Dimetilcarbonato+x 2 Hexano 

s(g E )' 

0.009 

0.008 

0.007 

0.005 

0.013 

0.006 


Ref.[15], Ref.[16], Ref.[17], 

s(Vi)' 

0.214 

0.192 

0.119 

0.105 

0.234 

0.106 


Ref.[18] 

s(/)" 

0.008 

0.008 

0.007 

0.007 

0.007 

0.007 



s(Vi)" 

0.064 

0.059 

0.069 

0.056 

0.082 

0.055 



s(h E ) 

121 

107 

113 

77 

76 

62 



s(v h ) 

37 

37 

34 

35 

34 

33 

8 

xiDimetilcarbonato+x 2 Octano 

s(/)' 

0.038 

0.034 

0.030 

0.022 

0.025 

0.022 


Ref.[15], Ref.[16], Ref.[17], 

s(Vi)' 

0.359 

0.336 

0.382 

0.205 

0.302 

0.177 


Ref.[18],Ref.[19] 

s(/)" 

0.027 

0.012 

0.028 

0.010 

0.031 

0.010 



s(Vi)" 

0.226 

0.265 

0.405 

0.057 

0.365 

0.039 



s(h E ) 

269 

111 

399 

80 

462 

85 



s(v h ) 

31 

30 

19 

20 

19 

20 



s(c£) 

2.1 

1.8 

1.9 

1.8 

1.8 

1.8 

9 

x\ Di meti lcarbonato+.r 2 Decano 

s(/)' 

0.058 

0.036 

0.056 

0.026 

0.056 

0.027 


Ref.[15], Ref.[16],Ref.[17], 

s(Vi)' 

0.702 

0.382 

0.718 

0.270 

0.632 

0.322 


Ref.[18],Ref.[20] 

s(/)" 

0.037 

0.031 

0.030 

0.032 

0.028 

0.032 



s(7ú" 

0.371 

0.222 

0.501 

0.162 

0.498 

0.135 
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s(h L ) 

s(v l ) 

s (c' v ) 

442 

62 

5.1 

109 

61 

5.1 

397 

18 

5.1 

84 

19 

5.1 

386 

14 

5.1 

110 

14 

5.0 

Éster + Alcohol 


ip 

*iEtanoato de Metilo+x 2 Metanol 

i'tó 1 )' 

0.006 

0.006 

0.006 

0.006 

0.006 

0.005 

1 

Ref.[21], Ref.[22], Ref.[23], 

Hy,)' 

0.039 

0.038 

0.037 

0.036 

0.036 

0.034 


Ref.[24] 

S(g')" 

0.007 

0.007 

0.007 

0.007 

0.007 

0.007 



s(yi)” 

0.033 

0.033 

0.034 

0.034 

0.034 

0.034 



s(h u ) 

8 

8 

7 

6 

6 

6 



s(v) 

2 

2 

2 

2 

2 

2 

11 

.X|Propanoato de Metilo+Jt 2 Etanol 

s( g F )’ 

0.020 

0.018 

0.021 

0.022 

0.021 

0.021 


Ref.[25], Ref.[26], Ref.[27], 

■s(j’i)' 

0.362 

0.341 

0.298 

0.270 

0.277 

0.219 


Ref.[28], Ref.[29] 

*(/)" 

0.011 

0.011 

0.011 

0.011 

0.011 

0.012 



s(?íT 

0.034 

0.032 

0.036 

0.034 

0.035 

0.032 



s(h') 

33 

33 

18 

26 

21 

19 



s(v E ) 

4 

4 

1 

1 

1 

1 

Aromático + Alcano 


12 

X\ Benceno+jc 2 Hexano 

s(g')' 

0.003 

0.003 

0.003 

0.003 

0.003 

0.004 


Ref.[30], Ref.[31], Ref.[32], 

s(y\)‘ 

0.013 

0.013 

0.013 

0.013 

0.012 

0.013 


Ref.[38], Ref.[39], Ref.[43], 

s(gT 

0.004 

0.004 

0.004 

0.004 

0.004 

0.004 


Ref.[44], Ref.[45], Ref.[49], 

s(yi)** 

0.015 

0.016 

0.016 

0.015 

0.016 

0.016 


Ref.[58] 

s(h E ) 

31 

31 

31 

31 

31 

31 



s(v h ) 

24 

22 

23 

22 

22 . 

22 



s(c E ) 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

13 

x ] Benceno+x 2 Heptano 

•Sfe 1 )' 

0.036 

0.036 

0.035 

0.036 

0.036 

0.036 


Ref.[33], Ref.[34], Ref.[35], 

s(y¡)‘ 

0.117 

0.117 

0.117 

0.118 

0.118 

0.118 


Ref.[36], Ref.[40], Ref.[41], 

s(g l )" 

0.016 

0.015 

0.016 

0.015 

0.015 

0.015 


Ref.[46], Ref.[50], Ref.[60] 

s(yi)** 

0.045 

0.043 

0.045 

0.043 

0.044 

0.044 



s(h E ) 

28 

31 

28 

28 

32 

28 



s(v h ) 

22 

22 

20 

20 

19 

19 



s(c;) 

0.1 

0.1 

0.1 

0.1 

0.0 

0.0 

14 

X\ Benceno+A: 2 Octano 

s(g k )' 

0.007 

0.020 

0.007 

0.016 

0.015 

0.015 


Ref.[31], Ref"[37], Ref.[40], 

s(y¡)' 

0.017 

0.040 

0.019 

0.031 

0.026 

0.026 


Ref.[42], Ref.[46], Ref.[50], 

s(gT 

0.014 

0.004 

0.017 

0.004 

0.008 

0.009 


Ref.[59] 

s(rT 

0.074 

0.048 

0.077 

0.058 

0.068 

0.069 



s(h E ) 

34 

29 

31 

24 

28 

28 



í(v E ) 

21 

20 

17 

17 

16 

17 



s( c; ) 

0.1 

0.1 

0.1 

0.0 

0.1 

0.1 

15 

x \ To 1 ueno+x 2 H exano 

s(g h Y 

0.011 

0.012 

0.012 

0.012 

0.011 

0.011 


Ref.[47], Ref.[48], Ref.[49], 

s(y\) 

0.024 

0.026 

0.023 

0.024 

0.024 

0.024 


Ref.[50], Ref.[61 ] 

s(gT 

0.020 

0.021 

0.021 

0.021 

0.020 

0.020 



s(y\)" 

0.054 

0.054 

0.054 

0.054 

0.053 

0.053 



sdñ 

19 

19 

18 

18 

17 

17 



i(v E ) 

8 

9 

5 

5 

6 

4 



s( c;) 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

16 

X]Tolueno+.r 2 Heptano 

s(TY 

0.012 

0.012 

0.011 

0.011 

0.011 

0.012 


Ref.[36], Ref.[50], Ref.[51], 

s(y,Y 

0.031 

0.030 

0.030 

0.030 

0.030 

0.030 


Ref.[52], Ref.[53], Ref.[54], 

s(gT 

0.005 

0.005 

0.005 

0.005 

0.005 

0.005 


Ref.[55], Ref.[56], Ref.[62] 

s(y¡)" 

0.030 

0.028 

0.029 

0.027 

0.029 

0.029 



s(h E ) 

12 

11 

11 

10 

10 

10 



■s(v') 

7 

7 

6 

7 

6 

7 



s(cp 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

17 

jCiTolueno+T 2 Octano 

s(g L Y 

0.011 

0.011 

0.011 

0.011 

0.011 

0.011 


Ref [48], Ref.[50], Ref.[57], 

s(y¡Y 

0.015 

0.014 

0.014 

0.013 

0.014 

0.014 


Ref.[58], Ref.[60] 

s(gT 

0.010 

0.010 

0.010 

0.011 

0.010 

0.011 



sqT 

0.019 

0.019 

0.018 

0.018 

0.018 

0.018 
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S(ff) 

12 

10 

12 

9 

13 

í(v E ) 

6 

6 

5 

5 

7 

s (cp 

0.0 

0.0 

0.0 

0.0 

0.0 


* Datos de £¿Fdeterminados en condiciones isotérmicas. ** Datos de üXFdeterminados en condiciones isobáricas 



Señalar que los valores de 5(g E ) y ¿'(^) presentados en la Figura 4 corresponden a 
valores medios de entre los obtenidos para el caso isobárico e isotérmico, diferenciados 
en la Tabla 2. Sin embargo, para el caso de las h E y de los v E se han reflejado los valores 
medios del coeficiente estadístico s 9 pero para estas propiedades obtenidos a varias tem- 
peraturas, cuando existen. 



Versión del Modelo 




Versión del Modelo 



Figura 4. Desviaciones estándar, s(y), con las 
diferentes versiones del modelo de la Tabla 1, 
aplicados en la correlación de los sistemas de 
la literatura. Los números indican la asigna- 
ción del sistema según la Tabla 2, y conside- 
rando el parámetro Ar¿j-libre. 
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En el segundo caso que se pretende analizar, se considera que el parámetro k X] 
depende de la presión y la temperatura de acuerdo a las propuestas realizadas para dicho 
parámetro en el Apéndice A.l. De esta forma, el modelo se acompaña de un mayor fun- 
damento físico al establecer en el procedimiento de correlación, las correspondientes 
funciones de dependencia para k\, k\¡ y k g , de acuerdo, respectivamente, a las ecuaciones 
(A.5), (A.7) y (A.ll). Para valorar la aplicación de los resultados obtenidos, estos se 
presentan de la misma forma que en la Tabla 3 y en la Figura 5. 

Ahora es preciso resaltar la excepción que se estableció en el procedimiento 
planteado para el pnmero, y se refiere a las correlaciones de c p E , por la falta de da- 

tos/funciones de c E (7). Aquí se consideró un valor constante del parámetro k\ } como 

cociente de las capacidades térmicas de los compuestos integrantes de la mezcla, 
k í .=c° 2 /cl l , similar forma de definición a la de los volúmenes para k y . 


Tabla 3. Error estándar de los ajustes conjuntos de datos de ELV en condiciones isotérmicas e isobáricas, de 
entalpías de exceso, h E i (J mor 1 ). volúmenes de exceso, lOV/Cm'-mol' 1 ). y capacidades térmicas de exceso. 
c p E /(J-mofK' 1 ) (cuando existen), mediante diferentes versiones del modelo (34) y calificados según la 

Tabla 1. En todos los casos se consideró al parámetro k¡¡ dependiente de la presión \ temperatura según las 
ecuaciones (A.5), (A.7) v (A.l 1) del apéndice, para evaluar respectivamente, los valores de k y . k h , y |- E . Para 
las correlaciones de se consideró un valor de k- x |= c p 2 / c pA . 

. 

Error 



Modelo 





estándar 

A 

B 

C 

D 

E 

F 

Alcohol + Alcano 


1 

Xi Etanol + x^ Hexano 

S(g\ 

0.032 

0.032 

0.021 

0.019 

0.020 

0.019 


Ref.[l], Ref.[2], Ref.[3], Ref.[4] 

s('A) 

2.611 

2.657 

1.664 

1.552 

1.255 

1.235 



s(gT 

0.029 

0.028 

0.020 

0.018 

0.017 

0.016 



s(y,)" 

1.691 

1.718 

1.167 

1.106 

1.077 

1.062 



s(h E ) 

158 

151 

123 

106 

109 

92 



sT) 

95 

95 

77 

80 

65 

67 

2 

xy Etanol - x-* Heptano 

sT) 

0.020 

0.020 

0.015 

0.014 

0.019 

0.018 


Ref.[l], Ref.[2], Ref.[3], Ref.[5] 

s(r,) 

0.401 

0.400 

0.282 

0.282 

0.296 

0.278 



s(gT 

0.021 

0.021 

0.013 

0.014 

0.011 

0.011 



s(y,)" 

2.002 

1.979 

1.594 

1.607 

1.428 

1.433 



s(h E ) 

179 

181 

151 

153 

136 

135 



s(v E ) 

98 

98 

76 

77 

56 

59 

3 

jr, Etanol - jr. Octano 

s(g E ) 

0.035 

0.033 

0.034 

0.033 

0.038 

0.037 


Ref.[l]. Ref.[2], Ref.[6] 

s(yd' 

1.161 

1.048 

0.899 

0.937 

0.S61 

0.867 



s(gT 

0.025 

0.023 

0.019 

0.019 

0.016 

0.015 



s(y-T 

0.488 

0.492 

0.391 

0.403 

0.354 

0.348 



s(h E ) 

213 

192 

186 

172 

156 

142 



s(\ f ) 

92 

91 

61 

62 

45 

51 

4 

ar¡ Etanol + r. N'onano 

s(gT 

0.028 

0.028 

0.024 

0.024 

0.022 

0.022 


Ref.[l] 

s(rT 

0.773 

0.776 

0.455 

0.445 

0.282 

0.378 



s(h E ) 

209 

212 

177 

178 

170 

169 



s(v E ) 

103 

102 

75 

75 

68 

65 

5 

Etanol + x** Decano 

s(gT 

0.040 

0.041 

0.030 

0.029 

0.029 

0.029 


Ref.[l] 

s(y¡)" 

0.467 

0.471 

0.250 

0.207 

0.206 

0.201 



s(h E ) 

228 

225 

186 

178 

177 

165 



s(v E ) 

132 

132 

113 

113 

96 

95 

6 

.r¡ Propan-1 -ol + r ; Heptano 


0.051 

0.035 

0.033 

0.025 

0.013 

0.012 


Ref.[7], Ref.[8], Ref.[9], Ref.[10] 

s(7\)' 

0.773 

0.551 

0.558 

0.399 

0.250 

0.222 



s(gT 

0.039 

0.028 

0.023 

0.020 

0.017 

0.014 



s(Y\)” 

0.411 

0.234 

0.369 

0.217 

0.326 

0.204 
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svñ 

i(v c ) 

199 

50 

196 

50 

174 

25 

168 

25 

151 

23 

147 

21 

Dialquilcarbonato + Alcano 


7 

x\ Dimetilcarbonato + x 2 Hexano 

5(g E )* 

0.016 

0.012 

0.012 

0.010 

0.012 

0.010 


Ref.[l 1], Ref.[12], Ref.[13], 


0.261 

0.208 

0.276 

0.236 

0.270 

0.234 


Ref.[14] 

S (/')’’ 

0.018 

0.019 

0.017 

0.017 

0.016 

0.017 



s(y ¡)” 

0.165 

0.155 

0.069 

0.082 

0.066 

0.082 



s(h E ) 

130 

117 

113 

91 

110 

80 



s(v h ) 

94 

94 

52 

57 

37 

38 

8 

X\ Dimetilcarbonato + X 2 Octano 

s(g E )' 

0.032 

0.033 

0.030 

0.031 

0.030 

0.031 


Ref.[l 1], Ref.[12], Ref.[13], 

s(yi)' 

0.351 

0.353 

0.353 

0.358 

0.344 

0.350 


Ref.[14],Ref.[15] 

s(gT 

0.011 

0.011 

0.010 

0.010 

0.010 

0.010 



s(y¡)" 

0.260 

0.262 

0.063 

0.063 

0.067 

0.070 



s(lñ 

148 

141 

117 

100 

115 

90 



s(v E ) 

122 

122 

82 

83 

53 

54 



s(c‘) 

6.961 

6.102 

6.828 

5.588 

6.799 

5.199 

9 

X\ Dimetilcarbonato +x 2 Decano 

s(g E r 

0.051 

0.051 

0.050 

0.046 

0.045 

0.046 


Ref.[ll],Ref.[12], Ref.[13], 

s(y¡)' 

0.673 

0.672 

0.475 

0.467 

0.439 

0.413 


Ref.[14],Ref. [16] 

S(g E )" 

0.032 

0.032 

0.029 

0.030 

0.029 

0.028 



s(y¡)" 

0.251 

0.245 

0.223 

0.215 

0.203 

0.192 



s(h E ) 

184 

174 

137 

115 

135 

108 



s(v E ) 

165 

165 

80 

104 

50 

62 



s(cp 

9.836 

8.504 

9.889 

8.156 

9.784 

7.908 

Éster + Alcohol 


10 

X\ Etanoato de Metilo +x 2 Metanol 

s(g E )‘ 

0.007 

0.008 

0.006 

0.006 

0.006 

0.005 


Ref.[17], Ref.[18], Ref.[l~9], 

s(y¡) 

0.058 

0.060 

0.052 

0.052 

0.049 

0.048 


Ref.[20] 

s (/)’’ 

0.005 

0.005 

0.004 

0.004 

0.004 

0.005 



■s(j’i)” 

0.034 

0.031 

0.037 

0.036 

0.037 

0.034 



s(h E ) 

50 

51 

44 

40 

36 

26 



s(v E ) 

8 

8 

3 

3 

3 

3 

11 

X[ Propanoato de Metilo +x 2 Etanol 

¡s(/)’ 

0.024 

0.022 

0.025 

0.023 

0.024 

0.023 


Ref.[21],Ref.[22], Ref.[23], 

í(’A)’ 

0.320 

0.305 

0.305 

0.275 

0.280 

0.261 


Ref.[24], Ref.[25] 

S(g h f 

0.012 

0.012 

0.013 

0.012 

0.012 

0.012 



¡s(}’¡)*‘ 

0.034 

0.030 

0.031 

0.033 

0.032 

0.033 



s(h E ) 

65 

44 

54 

28 

37 

21 



S(v £ ) 

8 

8 

2 

2 

2 

2 

Aromático + Alcano 


12 

X\ Benceno + x 2 Hexano 

s(g L )' 

0.006 

0.005 

0.006 

0.004 

0.006 

0.003 


Ref.[26], Ref.[27], Ref.[28], 

s(y¡)‘ 

0.032 

0.023 

0.029 

0.024 

0.029 

0.013 


Ref.[29], Ref.[35], Ref.[36], 

s(g E )" 

0.006 

0.004 

0.005 

0.004 

0.005 

0.004 


Ref.[40], Ref.[41], Ref.[42], 

s(y¡) 

0.036 

0.024 

0.029 

0.016 

0.029 

0.015 


Ref.[57] 

s(h E ) 

36 

34 

34 

31 

34 

31 



s(v E ) 

22 

22 

22 

22 

22 

22 



s(c E ) 

0.32 

0.10 

0.33 

0.06 

0.33 

0.05 

13 

x\ Benceno + x 2 Heptano 

s&y 

0.068 

0.054 

0.062 

0.043 

0.061 

0.039 


Ref.[30], Ref.[31 ], Ref.[32], 

s(y¡)' 

0.224 

0.145 

0.211 

0.139 

0.210 

0.142 


Ref.[33], Ref.[29], Ref.[37], 

s(g E )" 

0.034 

0.012 

0.025 

0.017 

0.025 

0.020 


Ref.[38], Ref.[43], Ref.[58] 

s(y¡)" 

0.137 

0.100 

0.096 

0.069 

0.088 

0.057 



s(h E ) 

38 

34 

32 

27 

32 

27 



s(v E ) 

27 

27 

20 

21 

20 

20 



s(cj¡) 

0.20 

0.09 

0.28 

0.02 

0.28 

0.03 

14 

jci Benceno + x 2 Octano 

s&Y 

0.040 

0.035 

0.038 

0.014 

0.038 

0.014 


Ref.[34], Ref.[27], Ref.[29], 

s(y¡)‘ 

0.118 

0.095 

0.115 

0.037 

0.117 

0.038 


Ref.[37], Ref.[39], Ref.[43], 

s(g E )" 

0.015 

0.007 

0.014 

0.021 

0.015 

0.019 


Ref.[57] 

■s(’/¡)” 

0.070 

0.043 

0.067 

0.074 

0.070 

0.070 



s(h E ) 

46 

43 

37 

29 

32 

25 
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v(/j' )/(J niol 1 ) - j f '¿ O O síx'/RT) 




J(V L ) 

s (c¡, ) 

27 

0.19 

27 

0.14 

17 

0.23 

16 

0.08 

16 

0.26 

17 

0.03 

15 

,T| Tolueno + Hexano 

s(¿‘)í 

0.041 

0.045 

0.044 

0.042 

0.042 

0.037 


Ref.[44], Ref.[45]. Ref.[46]. 

s(yif 

0.105 

0.097 

0.108 

0.095 

0.105 

0.087 


Ref.[47], Ref.[48], Ref.[59] 

s{g 1 )' 

0.028 

0.038 

0.034 

0.037 

0.032 

0.029 




0.101 

0.094 

0.091 

0.082 

0.087 

0.071 



s(h') 

19 

20 

20 

20 

20 

20 



i(v E ) 

10 

8 

7 

7 

6 

7 



s (c'„) 

0.05 

0.03 

0.05 

0.02 

0.05 

0.02 

16 

a*i Tolueno + x z Heptano 

s(g‘) 

0.018 

0.015 

0.015 

0.011 

0.014 

0.011 


Ref.[49], Ref.[50], Ref.[51], 

í( >’ i )'.. 

0.075 

0.051 

0.058 

0.033 

0.052 

0.031 


Ref.[52], Ref.[46], Ref.[53], 

s(g')‘ 

0.016 

0.010 

0.013 

0.006 

0.011 

0.005 


Ref.[33], Ref.[54], Ref.[48], 

s(y<j" 

0.052 

0.040 

0.036 

0.037 

0.029 

0.037 1 


Ref. [60] 

s(h l ) 

18 

17 

13 

12 

12 

11 



j(v E ) 

6 

6 

6 

6 

6 

6 



s (c{) 

0.21 

0.03 

0.23 

0.03 

0.23 

0.03 

17 

X\ Tolueno + x 2 Octano 

S(g'j 

0.011 

0.009 

0.010 

0.009 

0.011 

0.010 


Ref.[55], Ref.[46], Ref.[56], 

s(Y\) t 

0.024 

0.020 

0.015 

0.015 

0.015 

0.015 


Ref.[48]. Ref.[59] 

s(gT 

0.011 

0.013 

0.010 

0.011 

0.010 

0.011 



s(Y\)" 

0.024 

0.027 

0.021 

0.023 

0.020 

0.022 



S(h') 

14 

13 

12 

11 

12 

11 



J(V E ) 

9 

9 

8 

7 

7 

6 



s(c{) 

0.04 

0.04 

0.03 

0.03 

0.03 

0.02 


*Datos de £¿I 'determinados en condiciones isotérmicas. **Datos de £¿rdeterminados en condiciones isobáricas. 
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Figura 5. Desviaciones estándar, s(y), con las 
diferentes versiones del modelo de la Tabla 1, 
aplicados en la correlación de los sistemas de 
la literatura. Los números indican la asignación 
del sistema según la Tabla 3, y considerando el 
parámetro k^(p,T). 



0.0 


A 


B 


C 


D 


Versión del Modelo 


También, en la Tabla 3 se sombrearon las columnas correspondientes a los mo- 
delos que ocasionan las desviaciones mínimas globales de varias de las propiedades 
tratadas. 

La conclusión más inmediata es que la correlación de datos con el parámetro 
k\j(p 9 T) 9 como cabía esperar, produce peores resultados (cuantitativamente hablando) 
que cuando se considera Arij-libre, hecho que también se refleja al comparar las corres- 
pondientes gráficas, Figuras 4 y 5, de los parámetros de bondad de ajuste s, si bien las 
diferencias no se traducen claramente en las representaciones de las propiedades termo- 
dinámicas estudiadas. Las diferencias cuantitativas según las Tablas 2 y 3 no son sus- 
tancialmente relevantes. Por tanto, como resumen decir que al valorar el conjunto de 
los resultados de las diferentes versiones empleadas, tanto para este último caso como 
para el anterior, puede proponerse como aceptable, bajo un punto de vista práctico, el 
modelo calificado como “D” en la Tabla 1. 

Indicar por último que, con el fín de obtener una visión más cercana del trata- 
miento de datos de los sistemas presentados en las Tablas 2 y 3, en el Apéndice A.2 se 
realiza una aplicación numérica más detallada sobre el sistema binario (bence- 
no+hexano), el cual ha sido repetidamente utilizado como sistema estándar para la com- 
paración de las medidas experimentales realizadas con las diferentes propiedades y por 
diferentes autores. Además, coincide que es uno de los escasos sistemas de la literatura 
que ofrece un mayor número de propiedades determinadas experimentalmente y por 
tanto útil para comprobar la capacidad del modelo polinómico propuesto. 

7. Discusión de resultados 

En este trabajo se presenta un modelo matemático flexible, aplicable a la corre- 
lación de cualquier propiedad termodinámica de exceso, pero que se establece de mane- 
ra primaria para función de Gibbs de exceso de un sistema multicomponente, g E =g E (p, 

T, jci , X 2 , .); si bien, aquí se ha considerado su aplicación únicamente a sistemas bina- 

rios, que permite correlacionar simultáneamente datos de varias cantidades termodiná- 
micas. La flexibilidad del modelo da lugar a expresiones particulares para los casos iso- 
báricos o isotérmicos cuando interesen. 

La expresión genérica del modelo para un sistema binario resulta ser, 




( 47 ) 
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La aplicación del modelo a los sistemas encontrados en la literatura pennite llevar a 
cabo una elección del grado más adecuado del polinomio para los coeficientes g\(p,T), 
planteándose la expresión más extensa como: 

g, ( P.T )=g» +g¡z P'+gaPT- +^y+g i5 T 2 (48) 

De las relaciones formales termodinámicas, ver ecuaciones (1) y (2), se pueden 
obtener relaciones particulares para cada una de las propiedades de la solución. La for- 
ma polinómica del modelo (47) presenta la ventaja de seleccionar el grado del polino- 
mio acorde al número de cantidades termodinámicas que vayan a ser consideradas en el 
tratamiento matemático. 

Para correlacionar simultáneamente varias propiedades se empleó un procedi- 
miento de mínimos-cuadrados con un algoritmo genético, que minimiza las desviacio- 
nes de cada una de las propiedades que se consideren. De esta forma se determinan los 
valores óptimos para los coeficientes g \¡j de (48). Para el caso del parámetro k\] se anali- 
zaron dos posibilidades, una considerándolo como un parámetro adicional en el proceso 
de ajuste, que denominamos /r¡j-libre, y que ofreció en general las mejores correlaciones, 
y otro caso, en el que el parámetro tiene un valor determinado dado por las expresiones 
(A.5), (A.7) y (A.10), según se trate de correlacionar, respectivamente, las cantidades de 
v E , /? E y g E , y teniendo en cuenta la variación de dichas cantidades con la presión y la 
temperatura. Cuando existen datos de capacidades térmicas de exceso, c p E (7), el valor 

del parámetro se calcula como k c =Cp 2y /Cp, . 

Los resultados obtenidos indican que el modelo estima aceptablemente las dife- 
rentes propiedades de exceso para el conjunto de sistemas binarios escogidos para este 
estudio. Con referencia al modelo, si bien su forma polinómica contiene un elevado 
número de parámetros (un mínimo de quince con la versión “D”, pudiendo incluso au- 
mentar el número si los k\¡ se dejan libre en el procedimiento de ajuste), su expresión 
resulta la más adecuada, ya que por sucesivas derivaciones, que da lugar a su vez, a una 
subsiguiente reducción de parámetros, permite incluso aplicar un procedimiento de co- 
rrelación por etapas sucesivas en orden inverso a la derivación : Un ejemplo a este méto- 
do en etapas será (x,Cp )—>(x,/i L )—>(x,g E ), siempre y cuando existan los datos correspon- 

dientes a cada una de las propiedades. Además, permite una estimación aproximada de 
otras cantidades termodinámicas no obtenidas experimentalmente, como es el caso de la 
entropía, mejorándose la interpretación del comportamiento de la solución estudiada. 

En la práctica normal de la aplicación del modelo pueden despreciarse los térmi- 
nos correspondientes a la variación de la fracción activa “z” con la presión y la tempera- 
tura, ya que la contribución de los mismos supone un pequeño porcentaje del total, tal 
como se aprecia en el caso detallado en el apéndice, para el sistema binario estándar 
(benceno+hexano), con resultados aceptables en la correlación de todas las propiedades. 

Mencionar por último que el modelo propuesto en este trabajo es solo un primer 
paso para encontrar una solución más general, ya que la falta de una base de datos más 
extensa impide verificarlo en toda su extensión. Por ello, un densa labor de investiga- 
ción se presenta para un futuro a medio/largo plazo, ya que además de ampliar la base 
de datos existentes (que puede activarse con acuerdos entre varios grupos de investiga- 
ción), es necesario aún desarrollar otras aplicaciones del modelo, como su inclusión en 
la ecuación general de Gibbs-Duhem, su validación en la reproducción de diferentes 
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cantidades termofísicas como se recoge en las ecuaciones (1) y (2), y otras utilidades. 
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A. Apéndices 

A.l. Propuesta de cálculo para los parámetros k V] (p,T) 

Para complementar el estudio realizado es necesario conseguir una generaliza- 
ción de las ecuaciones particulares desarrolladas en trabajos anteriores para describir la 
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dependencia del parámetro ¿y con la temperatura [1-2] y con la presión [3], consideran- 
do en esta generalización la variación simultánea de una cantidad con la presión y la 
temperatura. Así, para un sustancia pura, si v(p,7), entonces: 


dv— 



'dv' 

— dT + 


UrJ p 

K d P, 


dp 


(A.l) 


dividiendo los diferentes términos por v y utilizando la definición establecida para los 
coefícientes de expansión isobárica a, y de compresibilidad isotérmica /?, se consigue la 
siguiente ecuación diferencial, 

dlnv-adT - pdp (A.2) 

cuya integración entre un estado inicial, tomado de referencia (po , Tq ), y otro cualquiera 
(p,7), da lugar a una ecuación que puede escribirse en función de densidades como si- 
gue: 


p(pj) = Po(Po’ T o)' e 


[p( P-Pq)-<*(T-T 0 )] 


(A.3) 


la cual se corresponde con una relación del tipo (f>(p,v,T) = 0, es decir, con una ecua- 
ción de estado en función de los coeficientes característicos de una sustancia pura. Aho- 
ra, la relación más general para k y (p,T), empleado en la correlación de volúmenes de 
exceso y establecido como cociente de volúmenes molares, tiene la forma siguiente: 


K(p,T)=k,(p,T)- 


V°(p,T)_M,p.(p,T)_ 

í Mj l 

r p-T 

V°(p,T) M,p,(p,T) 

l M \) 

K P'l’O > 


es decir: kj pj )=K{p n J^~^ n '^~ ai ' kT ' 


-ro)| 


JiP\-P'()p-P0)-{c\-<i'])T-T n )\ ^ 4 ) 


(A.5) 


estas relaciones permiten el cálculo del coeficiente K de acuerdo con la relación ya ob- 
tenida en anteriores trabajos [2,5] y que, en función de las dos variables, temperatura y 
presión, tiene la forma siguiente: 

K (P’ T )= ÍP’T )=K (p 0 ’T 0 (P’T F 73 (A.6) 


K (PI)=K (Po ’T 0 ÍKÍPo ’T 0 


(A.7) 


ya que los parámetros de área y ?\ no dependen de p y T. E1 uso de dicha ecuación 
está supeditada a conocer unos valores medios de a y /?para las sustancias puras en los 
adecuados intervalos de temperatura y presión. 

También, en los trabajos ya indicados antes se llevó a cabo la representación 
particular dela función de k x con la temperatura (caso isobárico) [1,5], y con la presión 
(caso isotérmico) [4]. Ahora se realiza la representación 3D con ambas variables y el 
resultado se recoge en la Figura A.l, donde se muestra una escasa dependencia de k y 
(para las mezclas consideradas) con la presión y la temperatura. 

De acuerdo a lo anterior, para el caso de la función de Gibbs de exceso, el pará- 
metro k g puede adoptar una forma análoga al utilizado para las otras cantidades termo- 
dinámicas analizadas. 0 sea, que podemos escribir una expresión del tipo: 
k g ( p,T )=k g (p 0 J 0 ]e^rTilp-Poy( a r a jl T - T o)\ (A.8) 

bajo un punto de vista práctico, al referimos al valor de k g (p 0 ,T 0 ) como una constante 


en unas condiciones de referencia, señalar que dichas condiciones no siempre coincidi- 
rán con aquellas para las que fueron obtenidos los coeficientes térmicos, por lo que la 
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ecuación (A.8) puede generalizarse aun más para escribirla como: 

k ^ p t) = d e ^ 1 -P') ^ p ~ p ° Lta ~ a i l T ~ T ° )J 

o bien: k e (p,T)=Ae [B(p - pohc{T - To)] 


(A.9) 

(A.10) 


0.25 



(a) 


’r 


(b) 



Figura A-l.- (a) Representación de la variación de los volúmenes molares V t ° , de compuestos puros 
(l:etanol, 2:hexano, 3:octano, 4:decano) con la presión y la temperatura. (b), Variación del parámetro 
de las mezclas binarias indicadas con los compuestos anteriores, con la presión y la temperatura. 

es decir, sustituyendo los coeficientes ténnicos, ay /?, y el valor de k g (p Q ,T 0 ) por unos 

parámetros a determinar en el procedimiento de ajuste, lo que originará un claro aumen- 
to del número de parámetros y posiblemente una ligera disminución de las desviaciones 
estándar, esto supone una sobreparametrización innecesaria en la mayor parte de los 
casos. 

Por ultimo, para el caso de las correlaciones de las capacidades ténnicas especí- 
ficas de exceso, c E , no se han encontrado valores de esta cantidad a diferentes tempera- 

turas y presiones, para los compuestos puros y para los sistemas binarios elegidos en 
este trabajo, por lo que se ha considerado un valor constante de /r¡j=Cp 2 /Cp, a una deter- 

minada presión y temperatura. 

A.2. Análisis de la función k l} (p^T) y su utilidad en el niodelo propuesto 

Como ejemplo de la utilidad del modelo propuesto, se presenta aquí una aplica- 
ción detallada de la versión D de la Tabla 1, ya que resultó ser bastante adecuada para la 
mayoría de los casos expuestos en la Tabla 2 y 3, indicando luego los resultados numé- 
ricos y las gráficas correspondientes al binario (benceno+hexano), por ser el sistema que 
presenta el mayor número de datos (isoténnicos e isobáricos) referenciados de varias de 
las propiedades, ELV, li , v E , y c E , utilizadas en el modelo. La ecuación completa de la 

versión del modelo “D” para la función de Gibbs de exceso, es la presentada en (34). 
Considerando las correspondientes definiciones establecidas con las ecuaciones (1) y (2) 
para conseguir las demás propiedades, pueden detallarse las relaciones particulares de 
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cada una de ellas. Como se comentó en la sección 5 de este artículo, en el tratamiento de 
datos termodinámicos de mezclas, es costumbre utilizar para el equilibrio entre fases la 
función de exceso de Gibbs en la forma adimensional, g E /RT , obtenida a partir de los 
coeficientes de actividad, y. Sin embargo, la primera derivada de esta fiinción con rela- 
ción a la temperatura para conseguir las entalpías de exceso, permite utilizar la forma 
no-adimensional expresada en (34). 

Así, la relación de g E /RT con la entalpía h E , viene dada por: 


h E 

~d{g E /RT)~ 


1 

Í5g E> | g E ’ 

g E T 

í 5gE l 

RT 2 

3T 

P,X 

RT 

^ 9T jp,x RT2 _ 

RT 2 RT 2 

IdT J 


(A.l 1) 

) p,x 

ecuación idéntica a la (36) o (38) que describe a la h E como función potencial de g E . Se 
trata ahora de definir cada uno de los dos sumandos de (A.ll), teniendo en cuenta la 
(A.7) para determinar ( dzjdT ), y considerar la expresión (37) para el término Y y su 
derivada. Con ello, la /? r , al considerar los coeficientes g\{p,T) de (34), resulta: 


/* E =z(l-z)¿ 
i=o V 


g\\+g\lp 2 + z \ L -gJ' 


dz 

~dT 


(A.12) 


Las ecuaciones para las otras cantidades termodinámicas de exceso, como y s E , pue- 
den conseguirse fácilmente desarrollando las respectivas ecuaciones. Así, la ecuación 


(40) para , se escribe ahora 


c I=-z{\-z) 2TY\ ^-+g\ 5 


dz 




/=o 


d{g,-gi- i)\/_7yf d2z ^ 


dT 2 


dT 

y teniendo en cuenta la (41), la entropía de exceso sería 

-s*=z{\-z)j^g,,p-^ + 2g, 5 T 

En el caso de utilizar datos de ELV isotérmicos, los volúmenes de exceso v E , y la fun- 
ción de Helmoltz a E , pueden determinarse, con las ecuaciones (43) y (44), respectiva- 
mente, dando lugar a: 


Iz'-tÍ— 

r \dT 


(A.13) 


(A.14) 


v E =z(l-z)^( 2 g¡2+g¡ 3 r)z i + 
i=o 

a E =z(l-z)¿ 
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dp 
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(A.16) 


Ahora es necesario un análisis sobre el empleo de los correspondientes paráme- 
tros Atíj y de sus derivadas con relación a p y T, ya que su intervención es clara en todas 
las ecuaciones anteriores. Resulta imprescindible estudiar su participación, y la del tér- 
mino donde se encuentre implicado, en el computo final de los resultados del modelo. 
Así, para el caso particular de la variación con la temperatura, es preciso considerar que:. 

(A.17) 
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pudiendo calcularse la variación del parámetro k con T de acuerdo a una de las expre- 
siones presentadas, (A.5), (A.7), (A.8), o (A.10). En caso de no considerar la variación 
de k con T, la nulidad del término (A.17) eliminaría los correspondientes en las ecuacio- 
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nes (A.12)-(A.14), con las correspondientes simplificaciones. La variación de z con p , 
ecuaciones (A.15), (A.16), se lleva a cabo con un procedimiento bastante análogo y no 
será expuesto para evitar repetición. 

A.2.1. Correlación de datos para sistema estándar (benceno+hexano) 

En esta sección se presentan los resultados de la aplicación detallada del modelo 
sobre el sistema (benceno+hexano), que utilizaremos para llevar a cabo un más comple- 
to análisis de resultados, su validación cuantitativa y cuantitativa, y verifícación si fuera 
necesario por otros autores. Como se ha indicado, una de las razones es que, para dicho 
sistema la bibliografia dispone de un buen número de datos experimentales para varias 
de las propiedades de mezcla aquí consideradas. Por ser sometido a una experimenta- 
ción continuada, los valores experimentales y/o correlaciones de las cantidades termo- 
dinámicas de este sistema son utilizados por muchos investigadores para su empleo co- 
mo referencia, o estándar de comparación, para la verificación de los procedimientos 
empleados. Concretamente, en la bibliografia consultada aparecen datos de ELV isobá- 
ricos [30] e isotérmicos [43-45], datos de h h [31,32,50] y de v r [38,39,50] a varias tem- 
peraturas, y datos de c p E [59]. 

Para llevar a cabo la correlación de datos se ha utilizado la versión del modelo 
indicado por “D” en la Tabla 1, y se ha seguido el procedimiento detallado en el esque- 
ma de la Figura 3. Se realizó la correlación simultánea de las cantidades cuyos datos 
experimentales para esa mezcla se recogen en la literatura indicada en el párrafo ante- 
rior y las ecuaciones (A. 12) a (A. 18). 

Otra apreciación importante sobre la forma polinómica del modelo, que podría 
dificultar a priori su aplicación, la constituye precisamente la consideración o no de la 
variación del parámetro /r¡j con la presión y la temperatura. Si se tuviese en cuenta esa 
variación, a través de ( dzldT ), aunque su inclusión sería científicamente más rigurosa, 
provocaría un aumento de los sumandos de cada una de las expresiones establecidas 
para las cantidades termodinámicas, ver ecuaciones (A.12) a (A.16), que de no ser así, 
quedarán reducidas sólo al primero de los sumandos.Por ello, en esta aplicación también 




Figura A-2. Para el sistema binario benceno+hexano. Representación de las contribuciones que 
aportan los diferentes sumandos separadamente, -parte superior, haciendo (dz/dT)=0; parte inferior, 
contribución del término con (dz/dT),- a la propiedad de exceso total: (a) sobre las h l según la 
ecuación (A.12); (b), sobre los v H según la ecuación (A.15). 
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hemos querido observar la contribución cuantitativa de esos términos en el computo 
total de la funciones termodinámicas derivadas de la g E . Así, gráficamente en la Figura 
(A-2) se aprecia que para el caso de las entalpías supone aproximadamente un 1,5% del 
total, mientras que para los volúmenes la participación es prácticamente nula. 

Por esta razón, los resultados de la aplicación sobre el sistema (benceno+ hexa- 
no) se presentan considerando solo el primero de los sumandos de cada una de las ecua- 
ciones (A.12) a (A.16). Los resultados finales se indican en la Tabla (A-l), distinguién- 
dose dos grandes columnas para la versión D elegida, dependiendo de la consideración 
realizada sobre el parámetro Ary. Así, el conjunto de la izquierda se refiere a los coefi- 
cientes g\(p,T) cuando el parámetro se considera libre en el procedimiento de ajuste, k\j- 
libre. Las tres últimas columnas de la derecha indican los resultados obtenidos de los 
coeficientes y las desviaciones al hacer variar los correspondientes parámetros k\¡ con la 
presión y temperatura, acorde a los datos experimentales utilizados; es decir utilizando 
valores de k(p,T) acorde a las relaciones obtenidas en el apéndice A.l. 

Se aprecia que las diferencias no son muy significativas en la diferente conside- 
ración del parámetro Ay, por lo que el empleo en este caso, de uno u otro camino en el 
proceso de correlación resulta indistinto. Por último, destacando una vez mas la buena 
correlación conseguida, al menos bajo un punto de vista cuantitativo, de las cantidades 
termodinámicas de este sistema binario, también deseamos hacer referencia a una repre- 
sentación en 3D y en 2D de los diferentes conjuntos de datos utilizados. En la Figura A- 
3 se muestran distintas gráficas con los puntos experimentales encontrados por diferen- 
tes autores y en distintas condiciones de presión y temperatura. También se presentan 
las curvas obtenidas en las correlaciones de las diferentes cantidades utilizando dos ver- 
siones del modelo para comparación, la versión tfc D” que hemos elegido aquí como base 
para mostrar el procedimento de trabajo, y con los parámetros de la Tabla (A-l) y A'y- 
libre, y la versión “F”, haciendo variar el parámetro k\ } con la p y T. Las diferencias en- 
tre ambos son inapreciables. 


Tabla A-l. Coeficientes, g\ } , para la versión del modelo “D” según Tabla 1, obtenidos en la correlación conjunta de y¡ 
de los ELV, isotérmico e isobárico, y las propiedades /t E , v E y c E del sistema binario *ibenceno+jc 2 hexano, (12) en 

las Tablas 2 y 3, junto a la desviación estándar s(M) de la cantidad M. Se considera el caso en que los parámetros k\- } - 
libre y el caso en que son funciones óepy T según las expresiones (A.5), (A.7), (A.10). 0corresponde ag^/RT. 


Parámetros k h , k g y k c determinados junto a los g¡j 

Parámetros k v , /: h , 

k c y k c función de p 

y T según las 


en el proceso de ajuste, es decir, Árj¡-libre. 


expresiones (A.5), (A.7), (A.10). 



gOj gl) 

g2} 

g0) 

SV) 

S2) 

Modelo D 

i = l 

-21.0 -889.2 

-782.8 

-1394.2 

833.4 

58.6 

j = 2 

2.267* 10‘ 8 8.355* 10' 9 

2.199><10' 9 

2.779x| 0' 8 

-I.l3|xl0' 10 

6.240x| 0' 10 

j = 3 

-8.608* 10' 9 -8.291* 10 9 

-1.605x10'' 

-1.140x10'* 

-4.396xl0‘ 9 

2.067x10' 10 

j = 4 

516069.1 146170.7 

209989.7 

914559.6 

-540801.0 

196343.6 

j = 5 

-3.917X10" 4 4.467* 10* 4 

2.393xl0' 3 

-2.084x10' 3 

7.449 xlO' 3 

-1.270xl0' 5 

s(0) 


0.004 UJ - 4 



0.004 u ' 3 ' 4 

s(X) 

k g = 1.124 

0.013 1 ' 2 ’ 34 



0.022 1 ' 2 ' 3 ' 4 

s(h E ) 

k h = 1.350 

3 1 5.6. 7 



3|5.6' 7 

•s(v E ) 

k v = 1.328 

227-8-9 



22 7 -8.9 

S(c E ) 

k c = 0.938 

0.08 10 



0.06 10 


'Ref. [30]; 2 Ref.[43]; 3 Ref.[44]; 4 Ref.[45]; 5 Ref.[31]; 6 Ref.[32]; 7 Ref.[50]; 8 R^[38]; 9 Ref.[39]; 10 Ref.[59] 
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Los resultados de la aplicación sobre el sistema (benceno+hexano) se presentan 
considerando únicamente el primero de los sumandos de cada una de las ecuaciones 
(A.12) a (A.16). Los resultados finales se indican en la Tabla (A-l), distinguiéndose dos 
grandes columnas para la versión del modelo elegida, dependiendo de la consideración 
realizada sobre el parámetro k\ j. Así, el conjunto de la izquierda se refiere a los coefi- 
cientes g\(p,T) cuando el parámetro se considera libre en el procedimiento de ajuste, /r¡¡- 
libre. Las tres últimas columnas de la derecha indican los resultados obtenidos de los 
coeficientes y las desviaciones al hacer variar los correspondientes parámetros k \j con la 
presión y temperatura, acorde a los datos experimentales utilizados; es decir utilizando 
valores de k(p,T) acorde a las relaciones obtenidas en el apéndice A. 1. 




Figura A-3. Representación de los valores experimentales de diferentes cantidades termodinámicas y las curvas de 
correlación obtenidas mediante la versión D, con los parámetros k g , k h , k x y k c determinados en el proceso de ajuste 

( - ), y la versión F del modelo con los parámetros k g , k hy k y y k c como función de p y T según las expresiones 

(A.5), (A.7), (A.10) (-) frente a la concentración del sistema binario XiBenceno+x 2 Hexano (12) de la Tabla 2 

y 3. (a) Función de energía de Gibbs de exceso del ELV isotérmico a T= 298.15 K (■), T = 328.15 K (Á) y T = 

333.15 K (•). (b) Función de energía de Gibbs de exceso del ELV isobárico a p = 101.32 ¿Pa (•). (c) Entalpía de 
exceso a T= 293.15 K (o), T= 298.15 K (A), T= 313.95 K (□), T =318.15 K (O) y T= 323.15 K (V). (d) Volumen 
de exceso a T= 288.15 K (o), T= 298.15 K (A), T= 303.15 K (□), T= 308.15 K (O), T= 318.15 K (V) y T = 

323.15 K(#). (e) Capacidad térmicade excesoa T= 298.15 K (o). 
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Se aprecia que las diferencias no son muy significativas en la diferente forma de 
considerar al parámetro por lo que utilizar uno u otro camino en el proceso de corre- 
lación resulta indistinto. Por último, destacar la buena correlación conseguida, al menos 
bajo un punto de vista cuantitativo, de las cantidades termodinámicas de este sistema 
binario, que se refleja en una representación en 3D y en 2D de los diferentes conjuntos 
de datos utilizados. En la Figura A-3 se muestran distintas gráficas con los puntos expe- 
rimentales encontrados por diferentes autores y en distintas condiciones de presión y 
temperatura. También se presentan las curvas obtenidas en las correlaciones de las dife- 
rentes cantidades utilizando dos versiones del modelo para comparación, la versión “D” 
que hemos elegido aquí como base para mostrar el procedimento de trabajo, y con los 
parámetros de la Tabla (A-l) y £jj-libre, y la versión U F”, haciendo variar el parámetro 
^ij con la /7 y T. Las diferencias entre ambos son inapreciables. 

Como resumen de esta aplicación, consideramos que el modelo utilizado en una 
correlación conjunta de las diferentes propiedades de una solución binaria resulta ade- 
cuado, ya que incluso permite una flexibilidad en la elección del número de parámetros 
acorde a la exigencia de los investigadores. Sin embargo, con la realización de este tra- 
bajo se ha mostrado la evidente necesidad de datos termodinámicos de soluciones en 
condiciones más extremas con el fín de contrastar aun más si cabe su utilidad. Por ello, 
la continuidad de este trabajo en un futuro permitirá ampliar la utilidad del modelo pro- 
puesto con la forma más óptima del mismo. 
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